Rotational Response of Two-Component Bose-Einstein Condensates in Ring Traps 
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We consider a two-component Bose-Einstein condensate (BEC) in a ring trap in a rotating frame, and show 
how to determine the response of such a configuration to being in a rotating frame, via accumulation of a 
Sagnac phase. This may be accomplished either through population oscillations, or the motion of spatial density 
fringes. We explicitly include the effect of interactions via a mean-field description, and study the fidelity of the 
dynamics relative to an ideal configuration. 

PACS numbers: 03.75.Mn, 06.30.Gv, 37.25.+k, 42.81.Pa 
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The Sagnac effect UJ is a rotational phenomemon de- 
scribing the phase shift, A^, between two coherent, counter- 
propagating waves traversing the same, closed path in a ro- 
tating frame. Originally discovered as an optical effect, it is 
actually more universal |2|; it has been observed in matter- 
wave interferometry experiments aiming to make precision 
measurements of rotation | 3 |, and has even been proposed as 
a method of testing general relativity |4 |. 

An attractive and theoretically simple geometry for observ- 
ing Sagnac-like effects in matter-waves is a ring trap, and 
with the excellent degree of control and precision now avail- 
able over magnetic and laser fields, the creation of such traps 
has recently been achieved by a number of groups worldwide 
(IH^I. Recent experiments with Bose-Einstein condensates 
(BECs) in ring traps 1 5 1 show how the coherent transfer of or- 
bital angular momentum to a trapped BEC O 23 can induce 
long-lived, superfluid flow. Two different flows (usually con- 
sidered to be counter-propagating, although this is not strictly 
necessary — as we will show, one of the flows may for ex- 
ample be zero) are required to observe the Sagnac effect in 
an atom-optical context. We show that there are a number 
of advantages in using a two-component BEC |11 1, made up 
of a single atomic species with two relevant internal states, 
particularly in ameliorating the frequently problematic effects 
of atom-atom interactions. It is therefore not necessary to 
assume negligible mean-field interactions fTT] to cleanly ob- 
serve the rotational response brought out by our proposed pro- 
tocols. We describe how the accumulation of a Sagnac phase 
can then be observed both through population oscillations be- 
tween the internal states, and by precession of density fringes 
within a particular internal state. We show how, in the case of 
density fringes, mean-field interactions from one component 
can stabilize the fringes in the other if the scattering lengths 
are approximately equal (as in, e.g., ^^Rb 1 13 1). Hence, within 
a mean-field picture, the repulsive interactions can in princi- 
ple be arbitrarily strong without affecting the interferometric 
signal. Finally, we discuss the sensitivity of our protocols. 

We consider a 2-component BEC composed of a single 
species with two relevant internal states, confined within an 
axially symmetric toroidal trapping potential considered to 
be insensitive to the internal state [Fig. [TJa)]. We employ a 
mean-field treatment, describing the sample with two coupled 
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FIG. 1: (color online), (a) Equipotential surface of an axially sym- 
metric toroidal trapping potential with radius p, in a frame rotating 
with angular frequency Q about the axis of symmetry (z axis), (b) 
In the quasi- ID limit, the equipotential surface shrinks to a ring, and 
the trapped atoms' spatial dynamics are in terms of the angle 6 only. 
Spatial fringes in the macroscopic atomic wavefunction will shift by 
AO = Q multiplied by the interrogation time. 



Gross-Pitaevskii equations: 
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where ^j(r) is the macroscopic wavefunction for atoms in in- 
ternal state j [normaHsed such that / <i^r|^Xr)P = 1], 
N is the total atom number, fioj is the energy difference be- 
tween the two internal states, the ajk are the ^-wave scatter- 
ing lengths, and M is the atomic mass. In terms of cylindri- 
cal coordinates {r, 6, z}, we take the confining potential to be 
y(r) = M[ajj(r - pf -h (jj\z^]I2, with torus radius p and an- 
gular trapping frequencies cOr, ojz- Assuming sufficiently tight 
radial and axial confinement the dynamics in these directions 
are "frozen out," permitting a quasi- ID description |[T4l [T5l 
[see Fig.[TJb)]. Projecting out the {r, z] dependences and mov- 
ing to a frame rotating counterclockwise about the z-axis with 
angular frequency transforms Eq. ([T]) to 
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where the time is now in units of r = Mp^ jfi, and frequencies 
are in units of r"^. The gjk = IMNpajk ^JoJrOJz/fi, and the 
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normalization condition is 2^=1 Jq^ = 1- 

To describe a BEC composed of atoms in a coherent inter- 
nal superposition state, we introduce the vector notation 
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We now present our Sagnac interferometry protocols. We 
assume all atoms to be initially in internal state 1 and in 
the motional ground state, such that the initial state i//^ has 
= 1/ V27r, = 0. Applying a resonant 7r/2 pulse (a "split- 
ting" pulse) to the internal two- state transition then yields 
n^^^^(0) = [/;,/2<A^ where 
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We imprint different angular momenta onto the spatial modes 
associated with the two internal states (e.g., by transfer of 
orbital angular momentum of light (91), producing ij/^^iO) = 
Uem^^l^iO), where 
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With this notation we can describe the symmetric case (m = 
0), the case where angular momentum is imprinted on only 
(m = £), and also permissable intermediate cases ({m, £} both 
either integer or half-integer). A free evolution f(T/2) follows 
[f(t) denotes an evolution governed by Eq. ^ for a time t, and 
T/2 is half the total interrogation time]. As the atom fields 
have uniform density, this takes a very simple form: 
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where ^1 = (m^-h^^)/2-h(gii -h2gi2 + g22)/87r, and ^2 = W2- 
2mi -h {g22 - gn)l^7i. We apply a n pulse Uj^, which swaps 
the two components, and allow another free evolution f(T/2) 
(completing the total interrogation time T). This negates the 
accumulated relative phase described by (fi2, producing 
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We repeat the angular-momentum-imprinting procedure U^m, 
which, due to the application of Uj^ dX t = r/2, undoes the 
relative diff'erence in angular momentum between the spatial 
modes associated with the two internal states. Following this 
by a second 7r/2 pulse (a "recombination" pulse) produces 
$\0) = U,i2UiJ^{e\ where 
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and there is no subsequent change to the popula- 
tions. We summarize this sequence by S n(T) = 
U,i2Ui^f{TI2)UJ{TI2)Ui^U,l2. The value of Q is in- 
ferred from the population, e.g., in internal state 2, N2 = 




FIG. 2: (color online), (a) State j mean relative population Nj/N 
as a function of interrogation time T, for (o) 7=1 and (x) j = 
2, following an S n(T) sequence, with ^ = 1. (b) Angular density 
fringes of state 1, with ^ = 1, m = 0, Q = 1, for (o) 7 = 0, (+) 
T = 0.1, ( ) r = 1, and (A) T = 8.3, following an S f(T) sequence. 
The interaction strengths gjk correspond to an ^^Rb configuration. 
All quantities are dimensionless. 



N[l - cos(2^Qr)]/2; the populations oscillate with period 
tt/^Q [see Fig. [2ja)]. Conveniently, any experimentally sig- 
nificant change of N2 from zero is a clear signal for finite Q. 
Note that the same response is obtained when angular momen- 
tum is imprinted on only (m = £) as for the symmetric case 
(m = 0). It is thus not essential to imprint angular momen- 
tum on both components, permitting a significant experimen- 
tal simplification. 

Alternatively, we may omit the second application of U^m, 
instead applying Uj:/2 directly to i//^(6) to produce 
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The atomic angular density in, e.g., internal state 1 is therefore 
N\il/i(e)\^ = A^{1 -h cos(2^[^ -h Qr])}/47r; the fringe spacing 
is 7i/£, and the fringe positions change with the total interro- 
gation time T with rate Q. Hence, measurable information 
about Q can be obtained without repeated angular momentum 
imprinting. If the measurement is not immediate, however, 
the if/j(0) continue to evolve. Simplifying to the case where 
,^11 =812 = 822 = for an initial condition 
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[e.g., formed from ij/^ by a t/;r/2f/^m^;r/2 sequence, or equiva- 
lent to Eq. ^ with 6 redefined and the global phase e'^"^^^ dis- 
carded] the mean-field contributions to Eq. ^ are g{co^^{i6)+ 
sin^(^^)]/27r = g/27r. A subsequent evolution f{t) yields 

^r(^) ^ ^-^ViV^-^(^+^0 / e^^-^comO + {^-m)n) \ ^^^^ 
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FIG. 3: (color online). Component 1 angular position density 
\il/i{e, t)\^ evolved by Eq. ^ from Eq. ([lo]) for m = 0, i^ = IQ = 0, 
gii = gii = 1000, with: (a) gn = 1000, (b) gn = 750, (c) gn = 500, 
(d) g\2 - 0. Note that |</^2(^)l^ is identical except for a n shift in 6, and 
finite Q causes the densities to precess by = Q.t. All quantities are 
dimensionless. 
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FIG. 4: (color online), (a) Fidelity F{t) = \ ^^"^ def{e,t)il/i{e,t)\^ of 
il/i{e, t) evolved by Eq. ^ from Eq. ([TO} with m = 0, ^ = 1, gn = 
g22 - 1000, for varying gu [(f>(6, t) is equivalent, but with gn = gn = 
gii]- (b), (c) Time-averaged fidelities F{T) = dtF{t), when 

gu = 1000, T = 10, for j = 1,2; ^^Rb lies within the (white) high- 
fidelity region. All quantities are dimensionless, and the maximum 
possible F(t) is 0.5. 



where the phases simplify to (fi = (n? -h £^)I2 -h glln, 
ip2 = oj/2 - 2ml. Hence, if the gjk are equal, the fringes in 
the two components stabilize each other, simply precessing 
around the ring with rate VI - m. Note also that the density 
fringes yielded by an S f(T) = f(T)Uj:/2U{mUj:/2 sequence 
are identical to those from a Uj:/2f(T/2)Uj:f(T/2)U£mUj:/2 se- 
quence when m = |[T6l [see Fig.[2jb)]. 

If the gjk are not equal, the fringe pattern can be strongly 
disrupted, as shown in Fig. [3] The scattering lengths can be 
very similar, however; in the case of ^"^Rb, if we consider F = 
2, rriF = I to be internal state 1, and F = I, mf = -1 to 
be internal state 2, then, in Bohr radii, an = 95.47, au = 
98.09, and ^22 = 100.44 (HI. Numerical results for an S f(T) 
sequence for various T, with gu = 974, gu = 1000, ^22 = 
1024, 1 = 1, m = 0, are shown in Fig.|2jb). The gjk correspond 
to an ^^Rb configuration with, e.g., p = 10"^ m, ojr = oJz = 
271 X 400 Hz, and N = 2800, which is consistent with current 
experimental capabilities |5], and for which the fringe profiles 
are only slightly perturbed ifTTll . 

Temporarily restricting ourselves to Q = 0, we note 
that Eq. ^ preserves initial periodic symmetry [i/^j(6) = 
i//j(6 + 27i/£)], and reflection symmetry [e.g., i//j(0) = 
(-1)^+ V/-^)]. and that, if if/jiO, t) form a solution to Eq. ([2]), 
then e~^^^^^^e^^^i//j(6 - mt, t) also form a solution. Further- 
more, if if/jiO, t) form a solution to Eq. (|2]) when II = 0, then 
if/jiO + Q.t, t) form a solution to Eq. (|2]) when Q 9^ 0. Hence, 
so long as the fringes remain resolvable [i.e., do not break up, 
as, e.g., in Fig. [3jd)], the fringe peaks precess with rate Q -m, 
and their form is independent of both Q and m. We may there- 
fore set D = 0, m = when considering the degree of fringe 



stabilization for unequal gjk numerically. In Fig.|4j for ^ = 1, 
we see that there is a broad region in parameter space exhibit- 
ing substantial stabilization, and Fig.[4jb) and Fig.Qc) show 
the ^^Rb parameter regime to be comfortably contained within 
this region. This result also indicates robustness to a compara- 
ble diflTerence in the local particle densities of the two internal 
states due to non-identical trapping potentials. Diflferent sta- 
bility regimes, and the rich dynamics shown in Fig. [3j could 
also be explored experimentally with the aid of a suitable Fes- 
hbach resonance |18|. 

As each atom is in a superposition state of the two inter- 
nal states, measuring the spatial distribution of atoms in, e.g., 
state 1 will on average project NI2 atoms into that internal 
state, with variance N/A. Hence, the standard deviation rel- 
ative to A/^ is 1 12 Va^, which can be considered negligible to 
the Gross-Pitaevskii level of approximation. An ideal in situ 
density measurement [6 1 specific to one internal state destroys 
the coherence between the internal states, but will in principle 
not aflTect the classical fields describing the positional states of 
the two components. The time-evolutions of the two compo- 
nents are still governed by Eq. ([2]), and so, evolving Eq. (10), 
if/iiO, t), 1/^2(0, t) continue to evolve according to Eq. \\\\ for 
the ideal gn = gu = g22 case 1 19 |. Consequently, assuming 
ideal, nondestructive measurements and gu ^ gu ~ ^22, the 
dynamics due to II 9^ may be tracked through repeated mea- 
surements within the same experimental run. Specifically, the 
fringe positions can be "zeroed" with a first measurement, and 
any subsequent precession monitored by later measurements. 
Finally, we note that, although we have assumed perfect ax- 
ial symmetry throughout, we do not expect the eflTect of any 
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potential asymmetries or corrugations to be significant if their 
scale is smaller than that set by the condensate chemical po- 
tential o. 

For an optical Sagnac interferometer, the fringe shift rela- 
tive to the fringe width Sl = AAQ^IAlc is commonly taken as 
a measure of the rotational sensitivity, where A is the interfer- 
ometer's enclosed area, is the optical wavelength, and c is 
the speed of light 1 1 1. The relevant angular shift of the spatial 
fringes we consider is = QJ, and the fringe widths scale 
as w = 7r/2^; hence d = AO/w = 21QT In. In the population- 
based protocol, d gives the number of instances the population 
alternates between and over a range [0, T] of interrogation 
times. It is instructive to now consider a more typical sim- 
ple Mach-Zehnder (MZ) configuration |20|. The de Broglie 
relations then yield for atoms of momentum p a wavelength 
A = h/p and velocity v = p/M, and hence Smz = 4AQM/h. 
Sensitivity of response therefore appears entirely determined 
by the enclosed area A, as opposed to the interrogation time 
T of our protocols. The time taken between the wavepacket 
splitting and its recombination in fact sets a natural timescale; 
in a ring geometry, this is given by Ti = np/v, which, us- 
ing V = M/pM and A = Trp^, becomes Ti = (AM/h)(27i/£). 
Hence, we may rephrase Smz = 2£QTil7i, which diff'ers from 
6 principally in that the time Ti at which it is possible to ex- 
tract useful information is fully determined by A and p, rather 
than being a free parameter [ 21 1. In an MZ configuration [201 
we expect v to be greater than that corresponding to small I in, 
e.g., ap = 10"^ m toroidal trap |5 1, implying moderate Ti for 
relatively large p. Our proposal is advantageous in using an 
intense, monochromatic source, where the usual associated is- 
sues of interatomic interactions have been circumvented, and 
where potentially useful information may be extracted at any 
time T. Large values of I [12J, although more challenging 
to generate, will also enhance the rotational response. Large 
p and N will aid in imaging; we note that ring traps of ra- 
dius p = 4.8 cm Q and ^^Rb condensates with ~ 10^ are 
achievable. The atomic shot noise also places a fundamental 
limit on the precision with which the spatial fringes or popu- 
lation oscillations can be measured, and hence Q inferred. 

In conclusion, we have considered the rotational dynam- 
ics of a single-species BEC, with two relevant internal states, 
within a ring trap configuration, and in a rotating frame. We 
have proposed a Sagnac-like interferometric protocol where 
the rotational sensing is manifest as time-dependent popu- 
lation oscillations, in a way that is insensitive to the atom- 
atom interactions arising within a mean-field picture. Simpler 
protocols involve observing the precession of density fringes 
around the ring. The fringes are robust for approximately 
equal interaction strengths (e.g., ^^Rb), and a range of strik- 
ing dynamics may also be observed by tuning the interactions 
with Feshbach resonances. All of these phenomena are ob- 
servable within the range of recent experimental advances. 
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